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Abstract 

It is suggested tliat gravity waves could, in several cases, be detected 
by means of already (or shortly to be) available technology, independently of 
current efforts of detection. The present is a follow-up on a recently suggested 
detection strategy based on gravity- wave- induced deviations of null geodesies. 
The new development is that a way was found to probe the waves close to the 
source, where they are several orders of magnitude larger than on the Earth. 
The effect translates into apparent shifts in stellar angular positions that 
could be as high as 10~^ arcsec, which is just about the present theoretical 
limit of detectability. 
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I. Introduction 



It was argued recently that gravitational waves could, in principal, be 
detected through a new non-mechanical effect [1,2]. This effect can be easily 
understood by replacing the usual textbook case of gravity-wave-induced 
variations in proper distances between initially static massive particles, by 
the analogous case for massless particles. One is then monitoring deviations 
in null, rather than in timelike geodesies. When the massless particles are 
photons, the timelike geodesies can be light rays coming from distant (e.g. 
stellar) sources, or perhaps from man-made space-based laser sources. Thus, 
what is to be measured in this approach, are gravity-wave-induced shifts in 
apparent angular separations of appropriate light sources. This effect should 
not be confused with the much smaller shifts in real angular separations of 
light sources, due to gravity waves actually hitting the light source [^. In 
0, it is the photons from the source (not the source itself) that encounter 
the waves, an that could happen arbitrarily close to the Earth. Nor should 
the effect be confused with yet another smaller shift, that due to the overall 
focusing of geodesies by the energy of the gravity waves. 

For an introduction to the physics of gravity waves and to the many 
proposed, and the few fully developed detection strategies, the reader is 
referred, in [4-22], to some reviews and key articles in the field. 

Consider, now, the following illustration of the effect, as discussed in 
1^: A gravitational pulse "with memory of position" [3,19,20,23-27] hits the 
Earth, coming from a supernova explosion, say, in the Northern sky. Alerted 
by the electromagnetic flash, the observer is now aware that a gravity wave 
has just whipped passed the Earth, and is traveling towards the Southern 
sky. For many years to come, light reaching the Earth from most sources in 
that hemisphere will have been through the gravity wave. Thus, if arbitrarily 
small angular shifts were detectable, the receding gravitational pulse could be 
"followed" for an extended duration. For most other schemes of gravity-wave 
detection, the same pulse would leave a single blip on the detector's record, 
and then would be lost for ever. Here, the gravitational pulse, all along its 
receding journey, keeps interacting with photons that eventually reach the 
observer. Hence, in a sense, the photons act as detectors that interact with 
the pulse farther and farther away from the Earth. 

At this point, the reader could be wondering what all the fuss is about, 
and she (or he) would be right, as far as this particular illustration is con- 
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cerned. For the waves get only dimmer (though at a neghgible rate,) as they 
keep moving away from us, and the expected apparent angular shifts (of the 
order of the wave's strength at the location of the photon-graviton interac- 
tion 0,) are over ten orders of magnitude smaller than what present day 
angular resolution astrometry can handle. 

However, as we shall now see, exploiting the same physics as above, but 
in a sense inverting the argument, one can arrive at a scheme for detecting 
gravity waves directly, using present-day resolution-power technology. 

Consider the following configuration. A steady source of gravity waves 
(e.g. a neutron star or a binary star, call it for short SG,) is giving off radi- 
ation with each polarization component having the form 

h= — exp{iQ{r - t)} , (1) 

where r is the distance to the center of the gravity-wave source SG, and Q 
is of the order of the time frequency of that source. (E.g. Q = 2uj when the 
source is a system with rotation frequency uj.) H is a constant that can be 
roughly estimated in the quadrupole approximation by 

H^MR^n"^ , (2) 

where M and R are the mass and size of SG. 

Now, consider a light source (dubbed SL hereafter,) situated even further 
away from us, along the line Earth-SG. (A situation where such a good 
alignment might happen naturally will be discussed below.) The photons we 
receive on the Earth from SL have, at some point in the past, crossed an area 
of relatively strong gravity waves. The maximum amplitude experienced by 
a photon is h^nax ~ H/b, where b is the impact parameter with respect to 
the center of SG. From the conclusions of 0, there is then some reason to 
suspect that the position of SL, as observed from the Earth, is continuously 
experiencing apparent angular shifts with an amplitude proportional to h{r = 
b). A priori, this effect could be larger than the one in by as many orders 
of magnitude as there is in the ratio r^j/fo, where ve is the distance from SG 
to the Earth. 

Upon inspection, it turns out that rather different considerations apply 
here than for the plane wave with memory of 0, leading to a modified 
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formula for the expected angular shifts. However, the calculations below 
show 1) that the effect here also is physical, 2) that its magnitude, in not too 
unlikely optimal situations, falls close to present limits of detectability. 

Before proceding to the next illustration, note that the angular shifts 
studied here are variations not about the actual position of SL, but about 
that position as displaced by non-radiative light deflections. However, those 
non-radiative deflections are either virtually constant in time (as in the ex- 
ample above,) and are therefore irrelevant here (besides being actually un- 
measurable,) or, as in the illustration below, they have a well-known time 
dependence and can be easily subtracted. In any case, Schwarzschild angu- 
lar shifts, for most interesting candidate sites, have a much larger time-scale 
than (and therefore should not interfere with potential measurements of) the 
expected shifts from gravity waves. 

Now, let our gravity-wave source SG be a rotating neutron star. A good 
candidate should be rotating rapidly enough to yield waves with a sufficient 
strength (see (2)), but not so rapid that the light crossing its gravity waves 
would shift too quickly for the sensitivity limit of the observer's hardware. 
(In the case of insufficient integration time, the overall effect of the waves 
would be a blurring of SL's image, with a caracteristic signature.) 

As for the role of the light source SL, let be played here by a companion 
star, locked with our neutron star SG into a binary system. Different SL 
candidates (ordinary main-sequence stars, pulsars, etc.,) will have different 
observational merits. Among the factors that come into play are the apparent 
magnitude of SL (shorter minimal integration times for larger magnitudes,) 
and its wavelength (better angular resolution for longer waves, best resolu- 
tion for radio waves.) Because quite different considerations apply for each 
type of companion SL, we limit ourselves here to deriving the main equations 
describing the effect, and obtaining orders of magnitude for reasonable opti- 
mal situations. A more detailed study of some promising actual astronomical 
sites, will be published elsewhere |^ . 

Taking a rapid pulsar as the companion SL (to the slower neutron star 
SG,) offers a particularly good opportunity for following the time evolution 
of the various parameters involved: Since it can be known to an excellent 
precision where each component of the binary pulsar is situated at any given 
moment (time resolution much finer than the period of either the neutron 
star SG or the pulsar SL,) one has an excellent time axis over which to draw 
eventual light deflection observations (see below.) 
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Just such a system has been considered recently to try and detect gravity 
waves through an effect that has httle to do with the angular shifts studied 



star's gravity waves on the pulsar's frequency that is considered. It is found 
to fall not too far from present limits of time-resolution power. In fact, 
the main problem there is less likely to be the time resolution than the 
photometric sensitivity, most pulsars being very dim radio sources (some 
notable exceptions are considered in [28].) 

Finally, we should perhaps call the attention on one more possible source 
of confusion. In the case (referred to two hereafter,) where the 

gravity-wave source SG is a neutron star and the light source SL an ordinary 
star or perhaps a pulsar, there are two quite different types of gravity waves 
involved. One is the neutron star SG's radiation, which has a frequency of 
typically several or many Hertz. The other is the gravitational radiation gen- 
erated by the binary system as a whole. These much longer waves (minimal 
period of a few hours,) are the ones probed in the case (case-one,) where the 
binary system as a whole was playing the role of SG, the light source SL 
being some much farther, unrelated body. In the following, when SG and SL 
are the two components of a binary (case-two,) the lower-frequency radiation 
of the whole binary system will not be considered. When we speak of the 
near-zone (distances less than one gravitational wavelength,) we mean the 
near-zone of SG, which is the neutron star in case-two, and the whole binary 
in case-one. 

This is not to say that it would not be interesting to consider the radiation 
of the whole binary even, in case- two. But, because we would then be not only 
in the near-zone of the source, but within the source itself, we would have 
to address the issue of disentangling dynamic Newtonian from relativistic 
contributions in that case [28]. 




In that paper, it is the Shapiro time-delay effect of the neutron 
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II. Equations for the effect 



Our first simplification will be to consider the effect of only one polariza- 
tion component of the gravity waves, namely the one that distorts geodesies 
in directions that are parallel to the plane of rotation of the gravity-wave 
source SG (that is the equatorial plane when SG is a neutron star, and the 
orbital plane when it is a binary.) This simplification does not affect substan- 
tially our conclusions, since the other component of the waves (the one with 
effects that are parallel to the rotation axis,) is usually of the same order of 
magnitude and frequency. (Although, attention to this other component can 
be extremely useful in disentangling radiative from dipole Newtonian effects 
in the near zone of SG.) Thus, we will be studying movements completely 
contained in the source's rotation plane. 

With these simplifications, the line element, in the spherical transverse- 
traceless gauge, has the simple form: 

df - dr^ - gdcl)'^ ^ , (3) 



where 

g = {l + hy . (4) 



Here the coordinates are centered on SG and light rays stretch from ~ 
to ~ vr. h is the gravitational- wave strength given by (1). As mentioned 
above, the Schwazschild component of the metric is not included, its con- 
tribution being either irrelevant (e.g. when SG and SL are very distant 
unrelated objects,) or easy to subtract from the observations (e.g when SG 
and SL are locked into a binary system.) For the fiat {h = 0) spacetime, 
of course, our equations should describe straight trajectories stretching from 
= to = TT, with closest encounter with SG at r = 6, where b is the 
"impact parameter". 

One can either stick to these polar coordinates or, alternatively, introduce 
a null coordinate v = r — t. It turns out that the gain from such a change 
for this particular problem is not considerable. So we shall stick here with t, 
r and (p. 
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Null geodesies are described by the equations 

dX 2^''^ 



(5) 



dp"" 1 







(6) 



-dX - 2'-P 








(7) 



Here = dx°'/dX are the photon's momenta. Commas indicate ordinary 
partial derivatives. 

Our aim is to determine as a function of r. We start by calculating 
p*^. We see immediately from (7) that, even for a finite h, there are straight 
line solutions, namely p'^ = 0. These are actually trajectories made of two 
(one incoming and one outcoming) components, with a discontinuous kink 
at r = 0. This is telling us that here, as in the plane wave case 0, photons 
that hit the wave fronts at right angles suffer no radiative deflection. 

For p^ 7^ 0, that is, for non-vanishing impact parameters, one gets from 
(7) the simple expression 



where L is an integration constant. 
Let us write (6) in the form 



where C is a constant to be determined shortly. Dividing this equation by 
(8) and using (1,4), one obtains 




(9) 



dr p"^ 
d(j) p'^ 



(1 + hy^ 



I 

JC 



■4> 2 + {l + iVLr)h 



d(f) 



(10) 



2(1 + h)r 
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In the absence of gravity waves, one should be able to recover the straight- 
line solution 



Tflat = -^—7 ■ (11) 
Sin0 

This fixes the value of C at 7r/2 , as one can see immediately from (10.11). 

Wc now switch to the radial variable u = 1/r. Then, deriving (10) with 
respect to 4>, and retaining only terms that are linear in h, we find 

u" + u = -h' ud(p - ^(u + in) , (12) 

Jn/2 2 

where primes indicate derivatives with respect to (f). 

In accordance with our linear treatment, the fiuctuation h can be ex- 
pressed in terms of (f> only by substituting, in (1), the expressions for r and 
t corresponding to the vacuum straight-line solution. These expressions are 
(11) and 

tflat = -0 -^—7 . (13) 

sm0 

(The origins of coordinates in (11,13) are chosen in such a way that tfiat — 
and Tflat = b at 4> = 7t/2.) Thus, we have 

u H ■ J. [ -or. 1 1 

/i = — sm exp oizo ; — > , (14) 

b [ sm0 J 

and 

/.'^^ = A^fcos0-zf)6l±^'| . (15) 
d(f) sm0 y sm0 y 

On the other hand, the integral in (12) is only needed to zeroth order. Hence, 
using (11), it is just — cosi 
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We are now in a position to write the equation governing light deflection 
from spherical gravity waves: 



u" + u 



6 sin I 



1 sin^ (J) — iQb 

2 



1 + cos sin ( 
sin (f) 2 



(16) 



This equation can be solved perturbatively by looking for a function ui 
satisfying 



U = Uo + Ui , 



(17) 



where Uq is the background solution, i.e., 



(18) 



which is the equation of a straight line (see (11)) . Then, the equation for ui 
is just (16), with u replaced by Ui on the left-hand side. 
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III. Interpretation 



The implications of (16) can be seen most clearly in the optimal case of 
very close encounter, when the impact parameter is less than one gravita- 
tional wavelength: b < A = 27r/fl. (Equation (16) can also be integrated 
analytically outside this region [32], but that will not be considered here.) 
As we shall see, such favorable configurations are not all that implausible 
(see also [28,32].) 

For what we called earlier case-one scenarii (SG and SL unrelated, far- 
apart objects,) if we consider a binary star like //-Sco (see below) as our 
gravity-wave source SG, then we have to find a candidate light source SL 
within an angle equal to SG's wavelength (about 18 light-hours) divided by 
SG's distance to the Earth (about 109 parsecs.) This gives an angular size of 
about 1 arcsec. This constraint is less stringent still if we consider a case-two 
situation (SG and SL tightly bound together.) Although the wavelength here 
is much shorter (say 0.1 light-second for a 10 Hz neutron star,) the distances 
are short enough (many such systems are only about 1 light-second across,) 
that the constraint on the orbit's inclination with respect to the Earth for 
obtaining a near- zone deflection is not too extreme (details below.) Among 
the over 500 known pulsars, several are known to be (and statistically many 
more should be) eclipsing or near to eclipsing binaries. In the logic of this 
paper, these obviously hold excellent prospects for gravity-wave detection. 

Going back to the inspection of (16) in the near- zone, a particular photon 
crosses the region of strongest h in less time than it takes the gravitational 
wave profile that he encounters to change appreciably. (Note that the 1/r 
fall-off of h makes it drop by a few orders of magnitude within a single 
wavelength). This means that the photon sees essentially a standing bulge 
of curvature (centered on the source), hence (16) can be integrated with 
h^H/r = i7 sin 0/6 (see (1,11)) . Subtracting this maximal trajectory from 
the one for when h is at its lowest gives the desired shift. 

Let us then compute the contribution of each term in (16) to Mi, and 
subtract the two values of at = (r = oo) to find the total shift for 
near-field defiection. We obtain 



H H TT TT 

A0 = ----M(- + l)+.(^- 
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+ J-J-^^H{^~^)+^^1 , (19) 

where we have kept the order of the terms in (16). 

First, we see that the effects of the f2- independent terms in (16) cancel 
out exactly. This means that there is no S c war scild- like light deflection. (In 
the Schwarzshild case, the RHS in (16) reduces to a term analogous to 
our hu oc sin^ 0.) 

Summing up all contributions, we obtain a formula for maximal shifts 
from near-field deflections: 



\A(PU,^^7rQH , (20) 



or 

3 2 



A^Ua^ ^ -IT hmax{r = A) . (21) 



Note that the amplitude does not increase as one penetrates deeper into the 
near-zone. This is because, a smaller value of b corresponding to a smaller 
angle of incidence (angle between a light ray and a given spherical wave- 
front), the effect (on A0) of the corresponding increase in h gets cancelled 
by the decrease in the sine of the angle (see [@.) 

As mentioned earlier, choosing a candidate site involves weighing the 
sometimes contradictory effects of a number of parameters, and finding a 
compromise that best accommodates the observational constraints. One ob- 
vious problem is that longer gravitational wavelengths, which correspond to 
wider near-zones, and hence to less stringent constraints on the alignment, 
often correspond to weaker gravity- wave amplitudes (see (2)). While we leave 
for 1^ a more detailed evaluation of various candidates' merits, we try here 
to estimate typical orders of magnitude for this effect in optimal, but realistic 
situations. 

Consider first, as our SG, the binary of giant stars /i-Sco in the Scorpius 
constellation. It is situated at a distance from the Earth te = 109 parsecs, 
and has an expected gravity wave frequency / ~ 1.6 x 10~^Hz. The expected 
dimensionless amplitude of gravity waves at r = r^; is h{rE) ~ 2 x 10"^*^, 
implying that H ^ 6cm. Such a source, which is among the strongest in 
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the Earth's vicinity, has unfortunately no chance of being detected by Earth- 
bound experiments, due chiefly to seismic noise. 

Let us see, in this instance, how far our effect falls from present limits 
of resolution power astrometry. The best resolution achieved today comes 
from Very Long Baseline Interferometry. There, the limit comes essentially 
from the shortest wavelengths that survive water and oxygen absorption (say 
0.1mm,) and the size of the Earth. The ratio gives about 10~^ arcsec. There 
are however definite projects for reaching 10~^ arcsec in a few years with 
space-based interferometry. 

Now, the numbers for /i-Sco imply a near-zone which spreads over an 
angular distance oi c/rf ~ 1 arcsec, where c is the speed of light. This 
does not seem too tiny to allow for a non-negligible chance of discovering 
a potential candidate for the role of SL. (Note that the Hipparcos sateUite 
routinely maps the sky at a resolution of 10~^ arcsec.) For obvious reasons, 
only the binaries that are closest to the Earth have so far receaved attention 
in the literature. But the closeness of SG to the Earth is not a sine qua non 
in the present context. Hence, a great many more binaries than the ones 
catalogued so far could be observationally interesting here. 

The shift's amplitude expected from (20) is 

\^(t>\max ~ ^tt'^/A 2 X lO'^arcsec , (22) 

which falls short of the limit of detectability by only one order of magnitude. 

The situation for two scenario (SG and SE in a binary) turns out to 

be even more promising. The period and strength of neutron-stars as gravity- 
wave sources are poorly known. The limits on neutron-star parameters are 
mostly drawn from pulsar data, but it is still uncertain how a typical ordi- 
nary neutron star compares to a typical pulsar. For instance, neutron-star 
rotation frequencies are thought to range all the way up to about lO^iJ^;, the 
observed higher limit for pulsars, but the lower limit is very uncertain.) Also, 
neutron-star gravity-wave strengths could vary widely depending on which 
of the possible mechanisms of gravity-wave production is actually at work. 
Eventually, only observation (why not of gravity waves,) will put stringent 
constraints on neutron-star parameters. 

To fix ideas, take SG to be a neutron star with / IHz, and SL to be a 
pulsar with a much higher frequency. SG and SL form a tight binary system. 
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Their average separation is only about one light-second. (The famous binary 
pulsar PSR 1913+16, e.g., is less than 3 light-seconds across, and several 
tighter binaries have been observed.) 

The alignment constraint here is really a constraint on the orbit's inclina- 
tion with respect to the Earth. For the hght from SL to cross the near-zone 
of SG on its way to us, the sine of that angle should be smaller than A/i?, 
where R is the orbit's half-size. But this ratio, given the numbers above, is 
typically close to one. Hence, we have here a set of candidates that naturally 
satisfy the alignment requirement. 

Now, if we assume that the neutron star is generating gravity waves 
through the CFS (Chandrasekhar- Friedman-Schutz) instability [30,31], then 
a possible value for H would he H ^ 10~^m (see [4] and references therein.) 
For / = IHz, this yields 



\^(t>\max ~ 10 '^arcsec . (23) 



And thus, even neutron-star gravity waves, some of the faintest to reach 
the Earth, could be detected afterall, and relatively soon. 

Finally, interesting candidates for the case-two scenario (SG and SL grav- 
itationally bound) are not limited to tight binaries of neutron stars. Take the 
X-ray source Her X-1. This is thought to consist of a neutron star (/ IHz) 
and a 2.2 solar-mass star, forming a binary system with an orbital period 
of 1.7 day. Because the orbit is larger here than in the case above (smaller 
A/i?), strong shifts such as (23) are more likely if the system is not too far 
from the eclipsing position. As it turns out, this is actually an eclipsing bi- 
nary, so excellent alignment occurs naturally every 1.7 day. We could also 
cite, in this class of candidates, several of the many multiple stars observed 
in the galaxy. These usually contain one or more binaries. The a-Gem sys- 
tem, for instance, is made of three binary stars, with periods of 0.8, 2.9 and 
9.2 days. Such configurations as this are, obviously, potential sites for rich 
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gravity-wave physics in the sense of this paper. 
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